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THE SLICE ALGORITHM FOR IRREDUCIBLE 
DECOMPOSITION OF MONOMIAL IDEALS 
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Abstract. Irreducible decomposition of monomial ideals has an increasing 
number of applications from biology to pure math. This paper presents the 
Slice Algorithm for computing irreducible decompositions, Alexander duals 
and socles of monomial ideals. The paper includes experiments showing good 
performance in practice. 
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1. Introduction 

The main contribution of this paper is the Slice Algorithm, which is an algo- 
rithm for the computation of the irreducible decomposition of monomial ideals. 
To irreducibly decompose an ideal is to write it as an irredundant intersection of 
irreducible ideals. 

Irreducible decomposition of monomial ideals has an increasing number of appli- 
cations from biology to pure math. Some examples of this are the Frobenius problem 
[U [2] , the integer programming gap [3] , the reverse engineering of biochemical net- 
works [4] , tropical convex hulls [5] , tropical cyclic polytopes [5] , secants of monomial 
ideals 6 , differential powers of monomial ideals [7] and joins of monomial ideals 
®- 

Irreducible decomposition of a monomial ideal / has two computationally equiv- 
alent guises. The first is as the Alexander dual of I [8], and indeed some of the 
references above are written exclusively in terms of Alexander duality rather than 
irreducible decomposition. The second is as the socle of the vector space R/I', 
where R is the polynomial ring that / belongs to and I 1 := I + (x\, . . . , x^) for 
some integer t >> 0. The socle is central to this paper, since what the Slice 
algorithm actually does is to compute a basis of the socle. 

Section [2] introduces some basic notions we will need throughout the paper and 
Section [3] describes an as-simple-as-possible version of the Slice Algorithm. Sec- 
tion 13.41 contains improvements to this basic version of the algorithm and Section 
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[5] discusses some heuristics that are inherent to the algorithm. Section [6] examines 
applications of irreducible decomposition, and it describes how the Slice Algorithm 
can use bounds to solve some optimization problems involving irreducible decom- 
position in less time than would be needed to actually compute the decomposition. 
Finally, Section [7] explores the practical aspects of the Slice Algorithm including 
benchmarks comparing it to other programs for irreducible decomposition. 

The Slice Algorithm was in part inspired by an algorithm for Hilbert-Poincare 
series due to Bigatti Et Al. [S]. The Slice Algorithm generalizes versions of the 
staircase-based algorithm due to Gao and Zhu [TU] (see Section 15. 2p and the Label 
Algorithm due to Roune [11] (see Section |5T5|) . 

2. Preliminaries 

This section briefly covers some notation and background on monomial ideals 
that is necessary to read the paper. We assume throughout the paper that /, J 
and S are monomial ideals in a polynomial ring R over some arbitrary field k and 
with variables x\, . . . ,x n where n > 2. We also assume that a, b, p, q and m are 
monomials in R. When presenting examples we use the variables x, y and z in 
place of x\ , X2 and X3 for increased readability. 

2.1. Basic Notions Prom Monomial Ideals. If v G N n then 

We define \fxf := x supp ^' where (supp(«))j := min(l, u-j). Define n (m) := 

such that e.g. 7T (x^ 1 ' 2 ' 3 ') = a* - ' 1 ' 2 ). 

The rest of this section is completely standard. A monomial ideal I is an ideal 
generated by monomials, and min (/) is the unique minimal set of monomial gen- 
erators. The ideal (M) is the ideal generated by the elements of the set M. The 
colon ideal I : p is defined as (m\mp G /). 

An ideal / is artinian if there exists ate N such that x\ G / for i — 1, . . . , n. 
A monomial of the form x\ is a pure power. A monomial ideal is irreducible if it 
is generated by pure powers. Thus ^x 2 ,y^) is irreducible while (x 2 yj is not. Note 
that (x) C y] is irreducible and not artinian. 

Every monomial ideal / can be written as an irredundant intersection of irre- 
ducible monomial ideals, and the set of ideals that appear in this intersection is 
uniquely given by /. This set is called the irreducible decomposition of /, and we 
denote it by irr (/). Thus irr ((a; 2 , xy, y 3 )) = {(x 2 , y) , (x, y 3 ) }. 

The radical of a monomial ideal / is \fl := (y'm|m£ min(/)). A monomial ideal 
I is square free if \fl — I. A monomial ideal is (strongly) generic if no two distinct 
elements of min (J) raise the same variable Xi to the same non-zero power [T2l I13j . 
Thus (x 2 y,xy 2 ) is generic while (xyz 2 ,xy 2 z) is not as both minimal generators 
raise x to the same power. In this paper we informally talk of a monomial ideal 
being more or less generic according to how many identical non-zero exponents 
there are in min (/). 

A standard monomial of I is a monomial that does not lie within /. The expo- 
nent vector v G N n of a monomial m is defined by m — x v — x\ x ■ ■ ■ x^ n . Define 
deg x (x v ) := Uj. We draw pictures of monomial ideals in 2 and 3 dimensions by 
indicating monomials by their exponent vector and drawing line segments separat- 
ing the standard monomials from the non-standard monomials. Thus Figure QJa) 
displays a picture of the monomial ideal (x 6 ,x 5 y 2 7 x 2 y A , y e ). 
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Figure 1. Examples of monomial ideals. 



2.2. Maximal Standard Monomials, Socles And Decompositions. In this 
section we look into socles and their relationship with irreducible decomposition. 
We also note the well known fact that the maximal standard monomials of I form 
a basis of the socle of R/I. 

Given the generators min (J) of a monomial ideal /, the Slice Algorithm computes 
the maximal standard monomials of /. We will need some notation for this. 

Definition 1 (Maximal standard monomial). A monomial to is a maximal standard 
monomial of I if to ^ I and mxi € / for i = 1, . . . ,n. The set of maximal standard 
monomials of / is denoted by msm (J) . 

The socle of R/I is the vector space of those to € R/I such that mxi — for 
i = 1, ... ,n. It is immediate that {to + I \m € msm (J) } is a basis of this socle. 

Example 2. Let / := (x 6 ,x 5 y 2 ,x 2 y 4: ,y 6 ') be the ideal in Figure QJ a). Then 
msm (J) = {x 5 y, x 4 y 3 , xy 5 } as indicated on Figure [2(b). Let J := (.x 5 y 2 , a; 2 y 4 ). 
Then msm (J) = {a; 4 y 3 } as indicated in Figure [2(c) . Finally, msm {(x 5 y 2 Yj = 0. 

We will briefly describe the standard technique for obtaining irr (/) from msm (/) 
[12]. Choose some integer t » and define 4>{x m ) = \rrii + 1 < i). 

Proposition 3 ([HI ex. 5.8]). The map is a bijection from msm (/ + (x\ , . . . , a^)) 
to irr (/). 

Example 4. Let / := (x 2 : xy) and /' := I + (a;*,?/) = (x 2 ,xy,y 3 ') where t = 3. 
Then msm(I') = {a;,?/ 2 } which </> maps to {(x 2 ,y) , (a;)} = irr (J). 

2.3. Labels. We will have frequent use for the notion of a label. 

Definition 5 (x^-label). Let d be a standard monomial of / and let to S min (I). 
Then to is an Xi-label of ci if m\dxi. 

Note that if to is an x^-label of d, then deg^fn) = deg^d) + 1. Also, a standard 
monomial d is maximal if and only if it has an a^-label rrii for i = 1, . . . , ri. So in 
that case dx± ■ ■ ■ x n = lcm" =1 to^. 

Example 6. Let / := (x 2 ,xz,y 2 7 yz,z 2 ') be the ideal in Figure [2](a). Then the 
maximal standard monomials of / are msm (/) = {xy, z}. We see that z has xz as 
an x-label, yz as a y-label and z 2 as a z-label. Also, xy has x 2 as an x-label and 
y 2 as a y-label, while it has both of xz and yz as z-labels. 
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Let J := I + (xy) be the ideal in Figure [2jt>) . Then msm (J) = {x,y, z}. Note 
that even though xy divides z ■ xyz, it is not a label of z, because it does not divide 
z ■ x, z • y or z • z. 




3. The Slice Algorithm 

In this section we describe a basic version of the Slice Algorithm. The Slice 
Algorithm computes the maximal standard monomials of a monomial ideal given 
the minimal generators of that ideal. 

A fundamental idea behind the Slice Algorithm is to consider certain subsets 
of msm (I) that are represented as slices. We will define the meaning of the term 
slice shortly. The algorithm starts out by considering a slice that represents all of 
msm (J). It then processes this slice by splitting it into two simpler slices. This 
process continues recursively until the slices are simple enough that it is easy to 
find any maximal standard monomials within them. 

From this description, there are a number of details that need to be explained. 
Section 13.11 covers what slices are and how to split them while Section 13.21 covers 
the base case. Section 13.31 proves that the algorithm terminates and Section 13.41 
contains a simple pseudo-code implementation of the algorithm. 

3.1. Slices And Splitting. In this section we explain what slices are and how to 
split them. We start off with the formal definition of a slice and its content. 

Definition 7 (Slice and content). A slice is a 3-tuple (I,S,q) where I and S 
are monomial ideals and q is a monomial. The content of a slice is defined by 
con (I, S, q) := (msm (I) \ S)q. 

Example [8] shows how this definition is used. 

Example 8. Let / := (a; 6 , x 5 y 2 , x 2 y A 1 y 6 ) and p :— xy 3 . Then / is the ideal de- 
picted in Figure [3](a) , where (p) is indicated by the dotted line and msm (I) = 
{x 5 y, x 4 y 3 , xy 5 } is indicated by the squares. We will compute msm (J) by perform- 
ing a step of the Slice Algorithm. 

Let Ii be the ideal / : p — (y 3 ,xy,x 4 ), as depicted in Figure E^b), where 
msm(Ii) = {x 3 ,y 2 } is indicated by the squares. As can be seen by comparing 
figures [3](a) and^b), the ideal I\ corresponds to the part of the ideal I that lies 
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Figure 3. Illustrations for example [8] 



within (p) . Thus it is reasonable to expect that msm ( I\ ) corresponds to the subset 
of msm (!) that lies within (p), which turns out to be true, since 

(1) msm (h)p — {x y 3 , xy 5 } = msm (I) n (p) 

It now only remains to compute msm (I) \ (p) . Let I2 ■= (x 6 , x 5 y 2 , j/ 6 ) as depicted 
in Figure [3Jc), where msm (I2) ■= {x 5 y,x 4 y 5 } is indicated by the squares. The 
dotted line indicates that we are ignoring everything inside (p). It happens to be 
that one of the minimal generators of J, namely x 2 y A : lies in the interior of (p), 
which allows us to ignore that minimal generator. We are looking at I2 because 

(2) msm (7 2 ) \ (p) = {x 5 y} = msm (I) \ (p) 

By combining Equation |T]) and Equation @ , we can compute msm (I) in terms of 
msm(Ji), msm (I2) and p. 

Using the language of slices, we have split the slice A := (I, (0) , 1) into the two 
slices A\ :— (Ii, (0) ,p) and A 2 :— (h, (p) ,1). By Equations |T]) and ((2|), we see 
that con (Ai) = msm (/) n (p) and con {A 2 ) = msm (J) \ (p). Thus 

con (A) = msm (J) = con (^4i) U con (A2) 

where the union is disjoint. 

Having defined slices and their content, we can now explain how to split a slice 
into two smaller slices. This is done by choosing some monomial p, called the pivot, 
and then to consider the following trivial equation. 



(3) 



(I, S, q) = ( con (I, 5, q) n (qp) ) U ( con (J, S, q) \ (qp) 



The idea is to express both parts of this disjoint union as the content of a slice. 
This is easy to do for the last part, since 

con (I, 5, q) \ {qp) = con (7, S+(p), q) 

Expressing the first part of the union as the content of a slice can be done using 
the following equation, which we will prove at the end of this section. 

msm (I) n (p) = msm (I : p) p 

which implies that (see Example [5]) 

con (I, S, q) n (qp) = con (J :p,S : p, qp) 



6 



BJARKE HAMMERSHOLT ROUNE 



Thus we can turn Equation ([3]) into the following. 

(4) con (7, S, q) = con (7 : p, S : p, qp) U con (7, S + (p) , q) 

Equation (|4]) is the basic engine of the Slice Algorithm. We will refer to it and 
its parts throughout the paper, and we need some terminology to facilitate this. 
The process of applying Equation ([4} is called a pivot split. We will abbreviate this 
to just split when doing so should not cause confusion. 

Equation (fj| mentions three slices, and we give each of them a name. We call 
the left hand slice (7, S, q) the current slice, since it is the slice we are currently 
splitting. We call the first right hand slice (I : p, S : p, qp) the inner slice, since 
its content is inside (qp), and we call the second right hand slice (I, S + (p) , q) the 
outer slice, since its content is outside (qp) . 

It is not immediately obvious why it is easier to compute the outer slice's content 
con (7, S + (p),q) than it is to compute the current slice's content con (7, S, q). The 
following equation shows how it can be easier. See Proposition 1111 for a proof. 

(5) msm (I)\S = msm (7') \ S, I' := (to G min (7) \tt (to) <£ S) 

This implies that con (7, S, q) = con (/', S, q). In other words, we can discard any 
element to of min (7) where ir (to) lies within S. We will apply Equation ([5]) when- 
ever it is of benefit to do so, which it is when ir (min (7)) n S ^ 0. This motivates 
the following definition. 

Definition 9 (Normal slice). A slice (I, S, q) is normal when ir (min (7)) PI S = 0. 

Example 10. Let I, p and I2 be as in Example[8] Then (7, (p) , 1) is the outer slice 
after a split on p. This slice is not normal, so we apply Equation [5] to get the slice 
(I2, (p) ,1), which is the slice A2 from Example See Figure [3] for illustrations. 

Proposition ! lll proves the equations in this section, and it establishes some results 
that we will need later. 

Proposition 11. Let J be a monomial ideal and let p be a monomial. Then 

(1) gcd(min(J)) divides gcd(msm (7)) 

(2) msm (I) n (p) — msm (I n (p)) 

(3) If p\ gcd(min (/)), then msm (I) = msm (I : p) p 

(4) msm (I) n (p) = msm (I : p)p 

(5) msm (7) \ S = msm (/') \ S, V := (to e min (I) |tt (to) g S) 

Proof. (pQ): Let d <G msm (I). Let U be an a^-label of d and let lj be an Xj-label 
of d where i ^ j. This is possible due to the assumption in Section [2] that n > 2. 
Then k\dxi and lj\dxj so gcd(min (7))| gcd(k, lj)\d. 
(|2|) : It follows from Lemma [12] below and (Q} that 

msm (J) n (p) = msm (In (p)) R (p) = msm (In (p)) 

If p\ gcd(min (7)) then p\ gcd(msm (7)) by |T]), whereby 

d G msm (7) 4=> {d/p)p I and (d/p)xip £ I for i = 1, . . . , n 

d/p ^ I : p and (d/p)xi € I : p for i = 1, . . . , n 

d/p £ msm (7 : p) <^4> d € msm (I : p)p 

(|4]): As p| gcd(min (7 n (p))) and (I C\ (p)) : p — I : p, we see that 

msm (7) n (p) — msm (7 n (p)) = msm ((7 f] (p)) : p) p — msm (I : p)p 
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([5]): Let d G msm (I)\S and let I G min (7) be an x^-label of d. Then I G min (/') 
since 7r (Z) \d 4- S. Thus dxi G I' since Z|dcci, so d G msm (/'). Also d £ I ~D I' . 

Suppose instead that d G msm(J') \ S. Then lix^ G 7' C 7. If <i G 7 then there 
would exist an m G min (!) \ min (/') such that m\d, which is a contradiction since 
then S 9 7T (m) |m|d ^ S*. Thus d £ I whereby d G msm (7). □ 

Lemma 12. Let A, 7? and C be monomial ideals. Then A C — B C implies 
that msm (A) n C = msm (7?) n C. 

Proof. Let d G msm (A) n C. We will prove that d G msm (B). 
d <£ B: If d G 5 then rfe BnC- AnC but rf^ 1 

dajj 6 S: Follows from dxi G A and d G C since then ctei £ AnC = BflC. □ 

3.2. The Base Case. In this section we present the base case for the Slice Algo- 
rithm. A slice (7, S, q) is a base case slice if 7 is square free or if x% ■ • ■ x n does not 
divide lem (min (7) ). Propositions [TBI and [HI show why base case slices are easy to 
handle. 

Proposition 13. If x% ■ ■ ■ x n does not divide lcm(min (/)), then msm (I) = 0. 

Proof. If msm (!) ^ then there exists some d G msm (I) . Let m G min (7) be an 
x^-label of d. Then Xi\m, so a?i|m| lcm(min (7)). □ 

Proposition 14. If 7 is square free and 7 ^ (x±, . . . , x n ) , then msm (7) =0. 

Proof. Let 7 be square free and let d G msm (7) . Let G min (7) be an a^-label 
of d for i = 1, . . . , n. Then d — it (lcm™ =1 m^) = 1 so rrii = Xi. □ 

3.3. Termination And Pivot Selection. In this section we show that some quite 
weak constraints on the choice of the pivot are sufficient to ensure termination. Thus 
we leave the door open for a variety of different pivot selection strategies, which is 
something we will have much more to say about in Section [5l 

We impose four conditions on the choice of the pivot p. These are presented 
below, and for each condition we explain why violating that condition would result 
in a split that there is no sense in carrying out. Note that the last two conditions 
are not necessary at this point to ensure termination, but they will become so after 
some of the improvements in Section [3.41 are applied. 

p £ S: If p G S, then the outer slice will be equal to the current slice. 
p y£ 1: If p = 1, then the inner slice will be equal to the current slice. 
p g I: See Section l4~4l and Equation (JT]) in particular. 
p\-K (lcm(min (/))): See Section 1431 and Equation ([8|) in particular. 

If a pivot satisfies these four conditions, then we say that it is valid. Proposition 
[15] shows that non-base case slices always admit valid pivots, and Proposition 1161 
states that selecting valid pivots ensures termination. 

Proposition 15. Let (I,S,q) be a normal slice for which no valid pivot exists. 
Then 7 is square free. 

Proof. Suppose 7 is not square free. Then there exists an Xi such that xf\m for 
some m G min (7), which implies that Xi ^ 7. Also, Xi (f. S since Xi\n (m) and 
(7, S, q) is normal. We conclude that IS 0j valid pivot. □ 

Proposition 16. Selecting valid pivots ensures termination. 
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Proof. Recall that the polynomial ring R is nocthcrian, so it does not contain an 
infinite sequence of strictly increasing ideals. We will use this to show that the 
algorithm terminates. Suppose we are splitting a non-base case slice A := (/, S, q) 
on a valid pivot where A\ is the inner slice and Ai is the outer slice. 

Let / and g be functions mapping slices to ideals, and define them by the ex- 
pressions f(I,S,q) := S and g(I,S,q) := (lcm(min (/))). Then the conditions on 
valid pivots and on non-base case slices imply that f{A) C f(A\), f(A) C /(A2), 
17(A) C g(Ai) and g(A) C 3(^.2). Also, if we let A be an arbitrary slice and we let 
A' be the corresponding normal slice, then f(A) C f(A') and g(A) C g(A'). 

Thus / and g never decrease, and one of them strictly increases on the outer 
slice while the other strictly increases on the inner slice. Thus there does not exist 
an infinite sequence of splits on valid pivots. □ 

3.4. Pseudo-code. This section contains a pseudo-code implementation of the 
Slice Algorithm. Note that the improvements in Section l3~4l are necessary to achieve 
good performance. 

The function selectPivot used below returns some valid pivot and can be im- 
plemented according to any of the pivot selection strategies presented in Section 
[5] A simple idea is to follow the proof of Proposition [15] and test each variable 
xi, . . . ,x n for whether it is a valid pivot. If none of those are valid pivots, then I' 
in the pseudo-code below is square free. 

Call the function con below with the parameters (J, (0) , 1) to obtain msm(7). 

function con(i, S, q) 

let I' := (to G min (I) |vr (m) ^ S) 

if X\ ■ ■ ■ x n does not divide lcm(min (/')) then return 
if I' is square free and I' 7^ (x\, . . . , x n ) then return 
if I' is square free and I 1 = (x±, . . . , x n ) then return {q} 

let p := selectPivot (/', S) 

return con (7' : p, S : p, qp) U con (/', S + (p) , q) 

4. Improvements To The Basic Algorithm 

This section contains a number of improvements to the basic version of the Slice 
Algorithm presented in Section [3l 

4.1. Monomial Lower Bounds On Slice Contents. Let ql be a monomial lower 
bound on the slice (I, S, q) in the sense that ql\d for all d S con (/. S, q). If we then 
perform a split on I, we can predict that the outer slice will be empty, whereby 
Equation ([4J specializes to Equation ((6]) below, which shows that we can get the 
effect of performing a split while only having to compute a single slice. 

(6) con (I, S, q) = con (J :l,S :l, ql) 

Proposition 1111 provides the simple monomial lower bound gcd(min(/)), while 
Proposition 1 1 71 provides a more sophisticated bound. 

Proposition 17. Let (7, S, q) be a slice and let 1(1) := lcm™ =1 where 

k ■= — gcd(min (I) C (x;)) 
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Then ql(I) is a monomial lower bound on (J, S, q). 

Proof. Let d £ msm ( I) and let m be an a^-label of d. Then Xi\m, so Uxi\m\dxi 
whereby k\d. Thus l(I)\d. □ 

Example 18. Let I :— (x 2 y,xy 2 ,yz, z 2 ). Then 1(1) = y and Equation ([6]) yields 

con (I, (0),l) = con(J:y, (0),y) 

where I : y = (x 2 , xy, z). As ^(7 : y) = x we can apply Equation ^ again to get 

con (7 : y, (0),y) = con((x,y,z), (0),xy) = {xy} 

We can improve on this bound using Lemma [20l below . 

Definition 19 (^-maximal). A monomial m 6 min(J) is Xi-maximal if 

< deg x .(m) = deg x .(lcm(min(J))) 

Lemma 20. Let d e msm (/) and let m be an Xi-\abe\ of d. Suppose that m is 
x j -maximal for some variable Xj. Then x% = Xj. 

Proof. Suppose that x% ^ Xj and let I be an Xj-\abel of d. Then 

deg x .(m) < deg Xj (d) < deg x .(l) < deg X3 (lcm(min (I))) = deg X] {m) □ 

Corollary 21. If m G min (/) is Xi-maximal for two distinct variables, then 
msm (I) = msm (/') where I' := (min (I) \ {m}). 

Corollary 22. Let (/, S, q) be a slice and let U :— ^- gcd(Mi) where 

Mi :— {m S min (/) \xi divides m and m is not a; j -maximal for any xj ^ Xi } 
Then q lcm™ =1 li is a monomial lower bound on (/, S, q). 

It is possible to compute a more exact lower bound by defining Muj\ and com- 
puting the gcd of pairs of minimal generators that could simultaneously be respec- 
tively Xi and jEj-labels. However, we expect the added precision to be little and 
the computational cost is high. If this is expanded from 2 to n variables, the lower 
bound is exact, but as costly to compute as the set msm(J) itself. 

Corollaries and fJU allow us to make a slice simpler without changing its 
content, and they can be iterated until a fixed point is reached. We call this 
process simplification, and a slice is fully simplified if it is a fixed point of the 
process. Proposition [23] is an example of how simplification extends the reach of 
the base case. 

Proposition 23. Let A := (I, S, q) be a fully simplified slice. If |min < n then 
A is a base case slice. 

Proof. Assume that x± ■ ■ ■ x n \ lcm(min (/)). Then for each variable Xi, there must 
be some rrii £ min (/) that is r^-maximal, and these to, are all distinct. Since 
I min < n this implies that min (I) — {mi, . . . , m n }. Thus Mi — {rrii} where Mi 
is defined in Corollary [22l Furthermore, since A is fully simplified, i gcd(Mi) = 1, 
so rrii — Xi and we are done. □ 

An argument much like that in the proof of Proposition 1231 shows that (/, S, q) 
is a base case if all elements of min (7) are maximal. If there is exactly one element 
m of min (I) that is not maximal, then one can construct a new base case for the 
algorithm by trying out the possibility of that generator being an Xi-label for each 
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Xi\m. One can do the same if there are k non-maximal elements for any k £ N, but 
the time complexity of this is exponential in A;, so it is slow for large k. 

Our implementation does this for k = 1,2, and implementing k — 2 did make 
our program a bit faster. We expect the effect of implementing k = 3 would be 
very small or even negative. 

4.2. Independence Splits. In this section we define /-independence and we show 
how this independence allows us to perform a more efficient kind of split. The 
content of this section was inspired by a similar technique for computing Hilbert- 
Poincare series that was first suggested in [TS] and described in more detail in [9] . 

Definition 24. Let A, B be non-empty disjoint sets such that AUB — {x±, . . . , x n }. 
Then A and B are I -independent if min (I) n (A) n (B) = 0. 

In other words, A and B are /-independent if no element of min (J) is divisible 
by both a variable in A and a variable in B. 

Example 25. Let I :— (x 4 , x 2 y 2 ,y 3 , z 2 , zt,t 2 y Then {x, y} and {z,t} are I- 
independent. It then turns out that we can compute msm (/) independently for 
{x, y} and {z, t}, which is reflected in the following equation. 

msm (I) = {x 3 yz : x 3 yt,xy 2 z,xy 2 t} = {x 3 y, xy 2 } ■ {z,t} 
= msm (/ n k[x, y\) ■ msm (I n n[z, i\) 

Proposition [26] generalizes the observation in Example l25l The process of apply- 
ing Proposition 1261 is called an independence split. 

Proposition 26. If A,B are /-independent, then 

msm (/) = msm (/ D k[A]) ■ msm (/ n k[_B]) 

Proof. Let A' := I n k[A] and B' := I H k[B\. If A' = (0) then msm (/) = by 
Proposition [T31 so we can assume that A' ^ (0) and B' ^ (0). It holds that 

min (/) = min (A 1 ) U min [B ) 

so for monomials a € k[A] and b £ k[B] we get that 

ab S / a G A' or b 6 S' 

and thereby 

ab^ I ^ A' and 6^5' 

which implies that 

a6 6 msm (/) 4^-ab ^ I and atei € / for Xi G A U /? 

A' and ax^ G A' for Xi £ A and 
6 ^ B' and € /?' for Xi £ B 
«a£ msm (A') and 6 £ msm (B') □ 

Given a slice (/, S, q), this brings up the problem of what to do about S when 
A and B are /-independent but not S'-independcnt. One solution is to remove the 
elements of min (S) n (A) n (-B) from min (S) when doing the independence split, 
and then afterwards to remove those computed maximal standard monomials that 
lie within (min (S) n (A) n (/?))■ Note that this problem does not appear if we use 
a pivot selection strategy that only selects pivots of the form x\ for t £ N. 
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Example 27. Let / be as in Example [25] and consider the slice (I, (x s y,y 2 z) , 1). 
Then y 2 z belongs to neither k[x, y] nor k[z, i], but we can do the independence split 
on the slice (I, (x 3 y) , 1) which has content {xy 2 } ■ {z 7 t} = {xy 2 z,xy 2 t} . We then 
remove (y 2 z) from this set, whereby con (J, (x 3 y,x 2 z), l) = {xy 2 t}. 

This idea can be improved by observing that when we know con (A', Sa, <?a)j 
we can easily get the monomial lower bound gcd(con (A 1 , Sa, 1a)), and we can 
exploit this using the technique from Section 14.11 This might decrease the size of 
min (S) n (A) n (B), which can help us compute con (B' , Sb,<1b)- 

This leaves the question of how to detect /-independence. This can be done in 
space OOn) and nearly in time 0(n |min (I)|) using the classical union-find algorithm 
[16] [17] I 1 ! See the pseudo-code below, where D represents a disjoint-set data struc- 
ture such that union(Z), Xi, Xj) merges the set containing x% with the set containing 
Xj. At the end D is the set of independent sets where D = {{xi, . . . , x n }} implies 
that there are no independent sets. The running time claimed above is achieved 
by using a suitable data structure for D along with an efficient implementation of 
union. See [HI [17] for details. 

let D := {{xi} , . . . ,{x n }}. 
for each m S min (7) do 

pick an arbitrary xi that divides m 
for each Xj that divides m do 
union(Z?, Xi, Xj) 

This is an improvement on the 0(n 2 |min algorithm for detecting indepen- 
dence suggested in [9]. That algorithm is similar to the one described here, the 
main difference being the choice of data structure. 

4.3. A Base Case Of Two Variables. When n = 2 there is a well known and 
more efficient way to compute msm (I). This is also useful when an independence 
split has reduced n down to two. 

Let {mi, . . . , TOfc} := min (I) where mi,. . . ,mk are sorted in ascending lexico- 
graphic order where x\ > X2- Let t(x u ,x v ) :— x^ Vl ' U2 '. Then 

msm (I) = {r(mi,m 2 ),r(m2,rn3), • . . ,r(mk-i,mk)} ■ 

4.4. Prune S. Depending on the selection strategy used, it is possible for the S 
in (/, S, q) to pick up a large number of minimal generators, which can slow things 
down. Thus there is a point to removing elements of min (S) when that is possible 
without changing the content of the slice. Equation ([7]) allows us to do this. 

(7) con(I,S,q)=con(I,S',q), S' := (m G min (S) \m £ I) 

Example 28. Consider the slice {(x 2 , y 2 , z 2 , yz) , (xyz) , 1). Then p := x is a valid 
pivot, yielding the inner slice ((x, y 2 , z 2 , yz) , {yz) , x). We can now apply Equation 
(0 to turn this into Ux, y 2 , z 2 , yz) , (0) , x). 

Proposition [TH] states that the Slice Algorithm terminates, and we need to prove 
that this is still true when we use Equation ([7]). Fortunately, the same proof can 
be used, except that the definition of the function / needs to be changed from 
/(/, S,q) — S to /(/, S, q) := I + S. Note that the condition on a valid pivot p that 
p (£ I is there to make this work. 

4t can also be done in space 0(n 2 ) and in time 0(n |min + n 2 ) by constructing a graph 
in a similar way and then finding connected components. 
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4.5. More Pruning of S. We can prune S using Equation {8j) below, and for 
certain splitting strategies this will even allow us to never add anything to S. 

(8) con (I, S, q) = con (/, S' , q) , S' := (m € min (S) \m divides tt (lcm(min (J))) ) 

To prove this, observe that any d € con (/, S, q) divides tt (lcm(min (/))). 

Example 29. Consider the slice ((x 2 , xy, y 2 ) , (0) , 1). Then p := x is a valid pivot, 
yielding the normalized outer slice ((xy, y 2 ) , (x) , 1). We can now apply Equation 
flU to turn this into ((xy, y 2 ) , (0) , 1). 

Similarly, Equation ([8]) will remove any generator of the form x\ from S. So if we 
use a pivot of the special form p = x\ , and we apply a normalization and Equation 
([H to the outer slice, we can turn Equation ((4]) into 

con (/, S, q) = con (I : p, S : p, qp) U con ((min (/) \ (pXi)), S, q) 

which for S = (0) and q — 1 specializes to 

msm (/) = msm (J : xf) x\ U msm ((min (/) \ 

An implementer who does not want to deal with S might prefer this equation to 
the more general Equation ((4|. 

We need to prove that the algorithm still terminates when using equations 
([7]) and |8]). We can use the same proof as in Proposition [TU except that we 
need to replace the definition of / from that proof with /(/, S, q) := I + S + 
(x" 1 , . . . ,x™ n ) where x u := lcm(min (/)). Note that the condition on a valid pivot 
p that p\tt (lcm(min (/))) is there to make this work. 

4.6. Minimizing The Inner Slice. A time-consuming step in the Slice Algorithm 
is to compute / : p for each inner slice (I : p, S : p, qp) . By minimizing, we mean the 
process of computing min (/ : p) from min (/) , which is done by removing the non- 
minimal elements of min (/) : p :— {m : p \m £ min (/) } where m : p := gcd ^ p ) ■ 

Proposition [30] below makes it possible to do this using fewer divisibility tests 
than would otherwise be required. As seen by Corollary I3T1 below, this generalizes 
both statements of [TBI Proposition 1] from p of the form x\ to general p0 See [HI 
Section 6] for an even earlier form of these ideas. 

Note that the techniques in this section also apply to computing intersections 
/ (p) of a momomial ideal with a principal ideal generated by a monomial, since 
min (/ (p)) — {lcm(m,p) \m G min (/ : p) }. 

The most straightforward way to minimize min (/) : p is to consider all pairs of 
distinct a, b G min (/) : p and then to remove b if a\b. It is well known that this can 
be improved by sorting min (/) : p according to some term order, in which case a 
pair only needs to be considered if the first term comes before the last. This halves 
the number of divisibility tests that need to be carried out. 

We can go further than this, however, because we know that min (J) is already 
minimized. Proposition 1301 shows how we can make use of this information. 

Proposition 30. Let x a , x b and x p be monomials such that x a does not divide x b . 
Then x a : x p does not divide x b : x p if it holds for i = 1, . . . , n that pi < V a,; < bi. 



2 This provides an answer to the statement from 1 181 p. 11] that "These remarks drastically 
reduce the number of divisibility tests, but they do not easily generalize for non-simple-power 
pivots, not even for power-products with only two indeterminates." 
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Proof. We prove the contrapositive statement, so suppose that x a does not divide 
x b and that x u := x a : x p divides x v :— x b : x p . Then there is an i such that asj > 
As max(pi, a^) = u,i + < + pi = max(j>i, 6,) we conclude that pi > a^. □ 

This allows us to draw some simple and useful conclusions. 

Corollary 31. Let a, b € min (7). Then a : p does not divide 6 : p if any one of the 
following two conditions is satisfied. 

(!) V® = V a : P 

(2) gcd(a,p)|gcd(6,p) 

Corollary 32. If a € min (J) and p\ir (a), then a : p £ min (7 : p) and a : p does 
not divide any other element of min (7) : p. 

We can push Proposition l30l further than this. Fix some monomial p, let x l := p, 
t 6 N", and define the function /, which maps monomials to vectors in N™, by 



(/(*"))< 



0, for i t = 0, 

I min(tij,£j + 1), for < s ; 7^ 0. 

Also, define the relation u -< v for vectors u, v € N n by 

u -< w if there exists an % € {1, . . . , n} such that ti > Ui > Vi. 

By Proposition [30l this implies that if a, b G min (7) and /(a) ^ then a : p 

does not divide b : p. This motivates us to define the function L(a) by 

L{a) := {m S min (7) |/(a) = /(to) } . 

If a, 6 € min (/) and f(a) -ft f(b), then no element of L(a) : p divides any element 
of L(b) : p. In particular, no element of L(a) : p divides any other. Note that the 
domain of L is a partition of min (I) . 

This technique works best when most of the non-empty sets L(a) contain con- 
siderably more than a single element, which is likely to be true e.g. if p is a small 
power of a single variable. Even in cases where most of the non-empty sets L{a) 
consist of only a few elements, it will likely still pay off to consider 7(1) and to 
make use of Corollary [521 



33. Let 7 : 


= (x 5 y,x 2 y 2 , 


x 2 z 3 ,xy 3 ,xyz 3 ,yz 2 } and 


p := x 3 . 


7(x 4 ) = 


{x 5 y} 


L(x 4 ):p=\ 


x 2 y} 


7(x 2 ) = 




L(x 2 ):p={ 




L{x) = 


{xy 3 ,xyz 3 } 


L(x):p={ 


y 3 ,yz 3 } 


7(1) = 


{y* 2 } 


L(l):p = \ 


yz 1 ) 



We will process these sets from the top down. The set L(x 4 ) is easy, since p\n (x 5 y) , 
so we do not have to do any divisibility tests for x 5 y. 

Then comes L(x 2 ). We have to test if any elements of L(x 2 ) : p divide any 
elements of L(x) : p or 7(1) : p. It turns out that x 2 y 2 : p\xy 3 : p and x 2 z 3 : 
p\xyz 3 : p, so we can remove all of 7(x) from consideration. We do not need to do 
anything more for 7(1), so we conclude that min (7 : p) = (x 2 y, y 2 , z 3 , yz 2 }. 
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4.7. Reduce The Size Of Exponents. Some applications require the irreducible 
decomposition of monomial ideals / where the exponents that appear in min (I) are 
very large. One example of this is the computation of Frobenius numbers [HH]. 

This presents the practical problem that these numbers are larger than can 
be natively represented on a modern computer. This necessitates the use of an 
arbitrary precision integer library, which imposes a hefty overhead in terms of time 
and space. One solution to this problem is to report an error if the exponents are 
too large, as indeed the programs Monos [TH] and Macaulay 2 [20] do for exponents 
larger than 2 15 — 1 and 2 31 — 1 respectively. 

In this section, we will briefly describe how to support arbitrarily large exponents 
without imposing any overhead except for a quick preprocessing step. The most 
time-consuming part of this preprocessing step is to sort the exponents. 

Let / be a function mapping monomials to monomials such that f(ab) — f(a)f(b) 
when gcd(a,6) = 1. Suppose that a\b =>■ f(a)\f(b) and that / is injective for each i 
when restricted to the set {x^ \x v S min (J) }. The reader may verify that then 

xi ■ ■ -a; n msm(/) = / _1 (:ei ■ • ■ x n msm ((/(min (J))))). 

The idea is to choose / such that the exponents in /(niin(J)) are as small as 
possible, which can be done by sorting the exponents that appear in min(7). If 
this is done individually for each variable, then |min (7)| is the largest integer that 
can appear as an exponent in /(min(/)). Thus we can compute msm(J) in terms 
of msm ((/(min (/)))), which does not require large integer computations. 

Example 34. If / := (x 100 , x 40 y 20 , y 9Q } then we can choose the function / such 
that (/(min(J))) = (x 2 ,xy,y 2 ). 

The underlying mathematical idea used here is that it is the order rather than 
the value of the exponents that matters. This idea can also be found e.g. in [T2l 
Remark 4.6], though in the present paper it is used to a different purpose. 

4.8. Label Splits. In this section we introduce label splits. These are based on 
some properties of labels which pivot splits do not make use of. 

Let (/, S, q) be the current slice, and assume that it is fully simplified and not a 
base case slice. The first step of a label split is then to choose some variable Xi such 
that min (I) n (x^ ^ {xi}. Let L := {x u £ min (/) \ui = 1}. Then L is non-empty 
since the current slice is fully simplified. Assume for now that \L\ = 1 and let I € L. 

Observe that if d G msm (J), then — \d if and only if / is an a^^-label of d, which 
is true if and only if Xi does not divide d. This and Equation (|3.1[) implies that 

con (J, 5,1) \ (xi) =con(/,5,l)n( — ) = con (l:—,S: — , — ) 

\ x% J y Xi x^ Xi j 

con (J, S, 1) n (x^ — con (/ : Xi,S : Xi,Xi) 

whereby 

con (I, S,q) = con (I : Xi, S : Xi,qXi) U con ( I : — ,S: — ,q — J 

\ Xi Xi Xi J 

This equation describes a label split on Xi in the case where \L\ = 1. In general \L\ 
can be larger than one, so let L = . . . , Ik} and define 



h _ / _ / h L_i \ \ I 
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for j = 1, . . . , k. Then con (Ij, Sj, qj) is the set of those d £ con (/, S, q) such that 
lj is an Xi-label of d, and such that none of the monomials ?i, . . . , lj-\ are ^labels 
of d. This implies that 

k 

(9) con (I, S, q) = con (I : x it S : x^qxi) [J con (Ij , Sj , qj ) 

i=i 

where the union is disjoint. This equation defines a label split on Xi. 

An advantage of label splits is that if / is artinian, S — (0) and \L\ = 1, then none 
of the slices on the right hand side of Equation ([9]) are empty. These conditions will 
remain true throughout the computation if the ideal is artinian and generic and we 
perform only label and independence splits. Example [35] shows that a label split 
can produce empty slices when \L\ > 1. 

Example 35. Let / := (x 4 , y A , z A , xy , x z). We perform a label split on x where 
?i := xy and I2 '■= xz, which yields the following equation. 

con (7, (0), 1) = con ((x 3 ,y, z), (0),x) (this is (I : x u S : x l ,qx l )) 

Ucon((x,y 3 ,z 4 ),(0),y) (this is (h,Sx,qx)) 

Ucon((x,y i ,z 3 ),(y),z) (this is (I 2 , £2, 92)) 

= {x 3 } U {y 3 z 3 } U = {x 3 ,y 3 z 3 } 

The reason that (I2, S2, 92) is empty is that both l\ and I2 are a;-labels of y 3 z 3 . 

Using only label splits according to the VarLabel strategy discussed in Section 
15.51 makes the Slice Algorithm behave as a version of the Label Algorithm [11] . See 
the External Corner Algorithm [2] for an earlier form of some of the ideas behind 
the Label Algorithm. 

5. Split Selection Strategies 

We have not specified how to select the pivot monomial when doing a pivot split, 
or when to use a label split and on what variable. The reason for this is that there 
are many possible ways to do it, and it is not clear which one is best. Indeed, it 
may be that one split selection strategy is far superior to everything else in one 
situation, while being far inferior in another. Thus we examine several different 
selection strategies in this section. 

We are in the fortunate situation that an algorithm for computing Hilbert- 
Poincare series has an analogous issue of choosing a pivot [9] . Thus we draw on the 
literature on that algorithm to get interesting pivot selection strategies [IB] , even 
though these strategies do have to be adapted to work with the Slice Algorithm. 
The independence and label strategies are the only ones among the strategies below 
that is not similar to a known strategy for the Hilbert-Poincare series algorithm. 

It is assumed in the discussion below that the current slice is fully simplified 
and not a base case slice. Note that all the strategies select valid pivots only. We 
examine the practical merit of these strategies in Section 17.21 

5.1. The Minimal Generator Strategy. We abbreviate this as MinGen. 

Selection. This strategy picks some element m € min (I) that is not square free 
and then selects the pivot n (m) . 
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Analysis. This strategy chooses a pivot that is maximal with respect to the property 
that it removes at least one minimal generator from the outer slice. This means 
that the inner slice is easy, while the outer slice is comparatively hard since we can 
be removing as little as a single minimal generator. 

5.2. The Pure Power Strategies. There are three pure power strategies. 

Selection. These strategies choose a variable Xj that maximizes |min (I) D (xi)\ pro- 
vided that xf | lcm(min (7)). Then they choose some positive integer e such that 
x1 +1 \ lcm(min(7)) and select the pivot x\. 

The strategy Minimum selects e := 1 and the strategy Maximum selects e := 
deg^. (lcm(min (/))) — 1. The strategy Median selects e as the median exponent of 
Xi from the set min (I) n (x$). 

Note that the Minimum strategy makes the Slice Algorithm behave as a version 
of the staircase-based algorithm due to Gao and Zhu [TU] , 

Analysis. The pure power strategies have the advantage that the minimization 
techniques described in Section 14.61 work especially well for pure power pivots. 
Maximum yields an easy inner slice and a hard outer slice, while Minimum does 
the opposite. Median achieves a balance between the two. 

5.3. The Random GCD Strategy. We abbreviate this as GCD. 

Selection. Let Xi be a variable that maximizes |min (!) n (xi)\ and pick three ran- 
dom monomials mi,m2,m,3 € min (!) n (xi). Then the pivot is chosen to be 
p := 7T (gcd(mi, ms, m^)). If p = 1, then the GCD strategy fails, and we might 
try again or use a different selection strategy. 

Analysis. We consider this strategy because a similar strategy has been found to 
work well for the Hilbert-Poincare series algorithm mentioned above. 

5.4. The Independence Strategy. We abbreviate this as Indep. 

Selection. The independence strategy picks two distinct variables Xi and Xj, and 
then selects the pivot p :— it (gcd(min (I) n (xiXj))). If p — 1, then the indepen- 
dence strategy fails, and we might try again or use a different selection strategy. 

Analysis. The pivot p is the maximal monomial that will make every minimal 
generator that is divisible by both Xi and Xj disappear from the outer slice. The 
idea behind this is to increase the chance that we can perform an independence 
split on the outer slice while having a significant impact on the inner slice as well. 

5.5. The Label Strategies. There are several label strategies. 

Selection. These strategies choose a variable Xi such that min (I) n (xj) ^ and 
then perform a label split on Xj. The strategy MaxLabel maximizes |min (/) n (xj)|, 
VarLabel minimizes i and MinLabel minimizes \ {x w e min (I) \vi = 1 }| while break- 
ing ties according to MaxLabel. 

Note that the VarLabel strategy makes the Slice Algorithm behave as a version 
of the Label Algorithm [TT] . 

Analysis. MaxLabel chooses the variable that will have the biggest impact, while 
MinLabel avoids considering as many empty slices by keeping |min [S)\ small. Min- 
Label is being considered due to its relation to the Label Algorithm. 
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6. Applications To Optimization 

Sometimes we compute a socle or an irreducible decomposition because we want 
to know some property of it rather than because we are interested in knowing the 
socle or decomposition itself. This kind of situation often has the form 

maximize v( J) subject to J € irr (J) 

where v is some function mapping irr (I) to M. We call such a problem an Irreducible 
Decomposition Program (IDP). As described in sections 16.31 and 1 6.41 applications 
of IDP include computing the integer programming gap, Frobenius numbers and 
the codimension of a monomial ideal. 

The Slice Algorithm can solve some IDPs in much less time than it would need 
to compute all of irr (I), and that is the subject of this section. Section luTTl explains 
the general principle of how to do this, while Section 16.21 provides some useful 
techniques for making use of the principle. Sections 16.31 and 16.41 present examples 
of how to apply these techniques. 

6.1. Branch And Bound Using The Slice Algorithm. In this section we ex- 
plain the general principle of solving IDPs using the Slice Algorithm. 

The first issue is that the Slice Algorithm is concerned with computing maximal 
standard monomials while IDPs are about irreducible decomposition. We deal with 
this by using the function 4> from Section 12.21 to reformulate an IDP of the form 

maximize v'(J) subject to J € irr (/') 

into the form 

maximize v(d) subject to d € msm (/) 

where v(d) := v'(4>(d)) and I := I' + (x\, ... , ) for some t » 0. 

It is a simple observation that there is no reason to compute all of msm (I) before 
beginning to pick out the element that yields the greatest value of v. We might as 
well not store msm (I), and only keep track of the greatest value of v found so far. 

We define a function b{I, S, q) that maps slices (I, S, q) to real numbers to be 
an upper bound if d £ con (J, S, q) implies that v(d) < b{I, S, q). We will now show 
how to use such an upper bound b to turn the Slice Algorithm into a branch and 
bound algorithm. 

Suppose that the Slice Algorithm is computing the content of a slice (I, S, q), and 
that b(I, S, q) is less than or equal to the greatest value of v found so far. Then we 
can skip the computation of con (J, S, g), since no element of con (7, S, q) improves 
upon the greatest value of v found so far. 

We can take this a step further by extending the idea of monomial lower bounds 
from Section I4TT1 The point there was that if we can predict that the outer slice of 
some pivot split will be empty, then we should perform that split and ignore the 
outer slice. That way we get the benefit of a split while only having to examine a 
single slice. In the same way, if we can predict that one slice of some pivot split 
will not be able to improve upon the best value found so far, we should perform the 
split and ignore the non-improving slice. The hard part is to come up with a way 
to find pivots where such a prediction can be made. Sections 16.31 and 16.41 provide 
examples of how this can be done. 

A prerequisite for applying the ideas in this section is to construct a bound b. It 
is not possible to say how to do this in general, since it depends on the particulars 
of the problem at hand, but Section I5~2l presents some ideas that can be helpful. 
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6.2. Monomial Bounds. In this section we present some ideas that can be useful 
when constructing upper bounds for IDPs of the form 

maximize v(d) subject to d G msm (7) . 

Suppose that v is decreasing in the sense that if a\b then v(a) > v(b). Then 
&(7, S, q) :— v(q) is an upper bound, since if d G con (7, S, q) then q\d so v(d) < v(q). 

Suppose instead that v is increasing in the sense that if a\b then v(a) < v(b). 
Then 6(7, S, q) := v(qir (lcm(min (/)))) is an upper bound, since if d G con (7, S, q) 
then d\q~K (lcm(min (/))) by Proposition 1361 below, so v(d) < u(g7r (lcm(min (/)))). 
Any monomial upper bound on con (7, S, q) yields an upper bound in the same way. 

Proposition 36. If d G msm (7) then d\ir (lcm(min (/))). 

Proof. Let d G msm (7) and let m, G min (7) be an x^-label of d for i = 1, . . . , n. 
Then d = tt (lcm™ =1 m,) divides n (lcm(min (7))). □ 

Sections 16.31 and 16.41 provide examples of how these ideas can be applied. 

6.3. Linear IDPs, Codimension And Frobenius Numbers. Let r G M. n and 

define the function v r (x u ) := u-r. Then we refer to IDPs of the form (TIT)]) as linear. 

(10) maximize v r (d) subject to d G msm (7) 

It is well known that the codimension of a monomial ideal I' equals the minimal 
number of generators of the ideals in irr(7'). The reader may verify that this 
is exactly the optimal value of the IDP (fT0|) if we let 7 := \TP + (x\, . . . ,x%) 
and r = (1, . . . , 1), noting the well known fact that the codimension of an ideal 
does not change by taking the radical. This implies that solving IDPs is NP-hard 
since computing codimensions of monomial ideals is NP-hard [151 Proposition 2.9]. 
Linear IDPs are also involved in the computation of Frobenius numbers [1] [2] . 

Let us return to the general situation of r and / being arbitrary. Our goal in this 
section is to solve IDPs of the form (TTU)) efficiently by constructing a bound. The 
techniques from Section RT21 do not immediately seem to apply, since v r need neither 
be increasing nor decreasing. To deal with this problem, we will momentarily 
restrict our attention to some special cases. 

Let a G K>o be a vector of n non-negative real numbers, and define v a (x u ) : — u-a. 
We will construct a bound for the IDP 

maximize v a (d) subject to d G msm (I) . 

This is now easy to do, since v a is increasing so that we can use the techniques 
from Section W% Specifically, v a (d) < v a (qTr (lcm(min (/)))) for all d G con (7, S, q). 

Similarly, let b G M< be a vector of n non-positive real numbers, and define 
Vb(x u ) :— u ■ b. We will construct a bound for the IDP 

maximize Vb(d) subject to d G msm (7) . 

This is also easy, since Vb is decreasing so that we can use the techniques from 
Section RT21 Specifically, Vb(d) < Vb(q) for all d G con (7, S, q). 

We now return to the issue of constructing a bound for the IDP (jTTJjl . Choose 
a G M™ and b G M™ such that r = a + b. Then we can combine the bounds for v a 
and Vb above to get a bound for v. So if d G con (7, S, q), then 

v(d) = v a (d) + v b {d) < v a (qir (lcm(min (7)))) + v b (q) =: b(I, S, q) 
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Now that we have a bound b, we follow the suggestion from Section 16.11 that we 
should devise a way to find pivots where we can predict that one of the slices will 
be non-improving. Let (I,S,q) be the current slice and let x u := lcm(min(/)). 

Suppose that fj is positive and consider the outer slice (I',S',q') from a pivot 
split on X4. We can predict that the exponent of Xi in our monomial upper bound 
will decrease from deg^ (q) + u» — 1 down to deg^. (q) . Thus we get that 

n(u l -l)<b(I,S,q)-b(l\S , ,q / ), 

whereby 

KJ',s',9')<^s.9)-nK-i), 

which implies that the outer slice is non-improving if 
(11) b(I,S,q)-ri(ui-l)<T, 

where r is the best value found so far. We can do a similar thing if is negative 
by considering the value of deg x (q') on the inner slice of a pivot split on . 

As we will see in Section ITU this turns out to make things considerably faster. 
One reason is that checking Equation (fTTjl for each variable Xi is very fast, because 
it only involves computations on the single monomial lcm(min (/)). Another reason 
is that we can iterate this idea, as moving to the inner or outer slice can reduce the 
bound, opening up the possibility for doing the same thing again. We can also apply 
the simplification techniques from Section 14.11 after each successful application of 
Equation (fTTj) . 

6.4. The Integer Programming Gap. Let c£ Q" and d € Z k , and let A be 

a k x n integer matrix. The integer 'programming gap of a bounded and feasible 
integer program of the form 

minimize c • x subject to Ax = d, x £ N" 

is the difference between its optimal value and the optimal value of its linear pro- 
gramming relaxation, which is defined as the linear program 

minimize c • x subject to Ax = d, x £ K™ . 

The paper [3] describes a way to compute the integer programming gap that in- 
volves the sub-step of computing an irreducible decomposition irr (/') of a monomial 
ideal Our goal in this section is to show that this sub-step can be reformulated 
as an IDP whose objective function v satisfies the property that a\b v(a) < v(b) 
whereby we can construct a bound using the technique from Section 16.21 

First choose t » and let I := /' + (a?* ) ■ • ■ j 2 -^ 1 ) so that we can consider 
msm (/) in place of irr (/'). Define tfj : N n i— > N n by the expression 

J Ui, for m < t, 
I 0, for Ui > t. 



So if t = 4 then ^(3, 4, 5) = (3, 0, 0). Define v(u) for u £ N n as the optimal value 
of the following linear program. We say that this linear program is associated to u. 

maximize c • (ip(u) — w) 

subject to -w) = 0, ioeI" 

and Wi > for those i where Ui < t 

The IDP that the algorithm from [3] needs to solve is then 

maximize v(u) subject to x u £ msm(/) . 
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By Proposition l37l bclow. we can construct a bound for this IDP using the technique 
from Section r6.2l Note that we can use this bound to search for non-improving outer 
slices for pivots of the form Xi in the exact same way as described for linear IDPs 
in Section [6751 

Proposition 37. The function v satisfies the condition that x a \x b =>■ v(a) < v(b). 

Proof. Let ej € N" be a vector of zeroes except that the i'th entry is 1. It suffices 
to prove that v(u) < v(u + ej) for u S N n . Let w £ K™ be some optimal solution to 
the linear program associated to u. We will construct a feasible solution w' to the 
linear program associated to u + ej that has the same value. We will ensure this by 
making w' satisfy the equation if>(u) — w = ip(u + ej) — w' . 
The case itj + 1 < t: Let w' := w + ej. 

The case itj + 1 = t: Let w' := w — Wjej. Note that the non-negativity con- 
straint on the i'th entry of w' is lifted due to Ui + 1 = t. 

The case itj + 1 > t: Let w' := w. Note that this case is not relevant to the 
computation since no upper bound will be divisible by □ 

7. Benchmarks 

We have implemented the Slice Algorithm in the software system Frobby [21], 
and in this section we use Frobby to explore the Slice Algorithm's practical perfor- 
mance. Section |7. II describes the test data we use, Section [7.21 compares a number 
of split selection strategies, Section 17.31 compares Frobby to other programs and 
finally Section 17.41 evaluates the impact of the bound optimization from Section [6] 

7.1. The Test Data. In this section we briefly describe the test data that we use 
for the benchmarks. Table [T] displays some standard information about each input. 
The data used is publicly available at |http : / / www . br oune . com/. 

Generation of random monomial ideals. The random monomial ideals referred to 
below were generated using the following algorithm, which depends on a parameter 
N e N. We start out with the zero ideal. A random monomial is then generated 
by pseudo-randomly generating each exponent within the range [0,7V]. Then this 
monomial is added as a minimal generator of the ideal if it does not dominate or 
divide any of the previously added minimal generators of the ideal. This process 
continues until the ideal has the desired number of minimal generators. The random 
number generator used was the standard C rand() function. 

Description of the input data. This list provides information on each test input. 

generic: These ideals are nearly generic due to choosing N = 30.000. 

nongeneric: These ideals are non-generic due to choosing N = 10. 

square free: These ideals are square free due to choosing N = 1. 

J51, J60: These ideals were generated using the reverse engineering algo- 
rithm of [3], and they were kindly provided by M. Paola Vera Licona. They 
have the special features of having many variables, being square free and 
having a small irreducible decomposition. 

smalldual: This ideal has been generated as the Alexander dual of a random 
monomial ideal with 20 minimal generators in 20 variables. Thus it has 
many minimal generators and a small decomposition. 
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1 
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40,920 


1 



Table 1. Information about the test data. 



t5d25p, t5d30p: These ideals are from the computation of cyclic tropical 
polytopes, and they have the special property of being generated by mono- 
mials of the form XiXj They were kindly provided by Josephine Yu. 

k4, k5: These ideals come with the program Monos [IH] by R. Alexander 
Milowski. They are involved in computing the integer programming gap of 
a matrix [3j. 

model4vars, model5vars: These ideals come from computations on alge- 
braic statistical models, and they were generated using the program 4ti2 
[2"2] with the help of Seth Sullivant. 

frobnl2dll, frobnl3dll: These ideals come from the computation of the 
Frobenius number of respectively 12 and 13 random 11-digit numbers pQ. 



7.2. Split Selection Strategies. In this section we evaluate the split selection 
strategies described in Section El Table [2] shows the results. 
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Table 2. Empirical comparison of split selection strategies. 



The most immediate conclusion that can be drawn from Table [2] is that label 
splits do well on ideals that are somewhat generic, while they fare less well on 
square free ideals when compared with pivot splits. It is a surprising contrast to 
this that the MinLabel strategy is best able to deal with J60. 

Table [2] also shows that the pivot strategies are very similar on square free ideals. 
This is not surprising, as the only valid pivots on such ideals have the form Xj, and 
the pivot strategies all pick the same variable. 

The final conclusion we will draw from Table [2] is that the Median strategy is 
the best split selection strategy on these ideals, so that is the strategy we will use 
in the rest of this section. The Minimum strategy is a very close second. 

7.3. Empirical Comparison To Other Programs. In this section we compare 
our implementation in Frobby [21] of the Slice Algorithm to other programs that 
compute irreducible decompositions. There are two well known fast algorithms for 
computing irreducible decompositions of monomial ideals. 

Alexander Dual [3 [23]: This algorithm uses Alexander duality and inter- 
section of ideals. Its advantage is speed on highly non-generic ideals. 
Scarf Complex [T2l [23] : This algorithm enumerates the facets of the Scarf 
complex by walking from one facet to adjacent ones. The advantage of 
the algorithm is speed for generic ideals, while the drawback is that highly 
non-generic ideals lead to high memory consumption and bad performance. 
This is because the algorithm internally transforms the input ideal into a 
corresponding generic ideal that can have a much larger decomposition. 
We have benchmarked the following three programs. 

Macaulay 2 version 1.0 [20]: Macaulay 2 incorporates an implementation 
of the Alexander Dual Algorithm. The time consuming parts of the algo- 
rithm are written in CH — h 
Monos version 1.0 RC 2 [19] : This Java program^ incorporates Alexander 
Milowski's implementation of both the Alexander Dual Algorithm and the 
Scarf Complex Algorithm. 
Frobby version 0.6 [21] : This C++ program is our implementation of the 
Slice Algorithm. 



There are two different versions of Monos that have both been released as version 1.0. We 
are using the newest version, which is the version 1.0 RC2 that was released in 2007. 
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How these programs compare depend on what kind of input is used, so we use 
all the inputs described in Section 17.11 to get a complete picture. In order to run 
these benchmarks in a reasonable amount of time, we have allowed each program 
to run for one hour on each input and no longer. Each program has been allowed to 
use 512 MB of RAM and no more, not including the space used by other programs. 
We use the abbreviation OOT for "out of time" , OOM for "out of memory" and 
RE for "runtime error" . 

The benchmarks have all been run on the same Linux machine with a 2.4 GHz 
Intel Celeron CPU. The reported time is the user time as measured by the Unix 
command line utility "time" . 

All of the data can be seen on Table El The data shows that Frobby is faster than 
the other programs on all inputs except for smalldual. This is because the Alexander 
Dual Algorithm does very well on this kind of input, due to the decomposition being 
very small compared to the number of minimal generators. The decompositions of 
J51 and J60 are also small compared to the number of minimal generators, though 
from the data not small enough to make the Alexander Dual Algorithm win out. 

It is clear from Table [3] that Macaulay 2 has the fastest implementation of the 
Alexander Dual Algorithm when it does not run out of memory. As expected, the 
Scarf Complex Algorithm beats the Alexander Dual Algorithm on generic ideals, 
while the positions are reversed on square free ideals. 

As can be seen from Table [3l the other programs frequently run out of memory. 
In the case of Macaulay 2, this is clearly in large part due to some implementation 
issue. However, the issue of consuming large amounts of memory is fundamental to 
both the Alexander Dual Algorithm and the Scarf Complex Algorithm, since it is 
necessary for them to keep the entire decomposition in memory, and these decom- 
positions can be very large - see frobnl3dll as an example. The Slice Algorithm 
does not have this issue. 

An advantage of the Slice Algorithm is that the inner and outer slices of a 
pivot split can be computed in parallel, making it simple to make use of multiple 
processors. The Scarf Complex Algorithm is similarly easy to parallelize, while 
the Alexander Dual Algorithm is not as amenable to a parallel implementation. 
Although Frobby, Macaulay 2 and Monos can make use of no more more than 
a single processor, multicorc systems are fast becoming ubiquitous. Algorithmic 
research and implementations must adapt or risk wasting almost all of the available 
processing power on a typical system. E.g. a non-parallel implementation on an 
eight-way system will use only 13% of the available processing power. 



7.4. The Bound Technique. In this section we examine the impact of using the 
bound technique from Section [6] to compute Frobenius numbers. 

Table |4] displays the time needed to solve a Frobenius problem IDP with and 
without using the bound technique for some split selection strategies. We have 
included a new selection strategy Frob that works as Median, except that it selects 
the variable that maximizes the increase of the lower bound value on the inner slice. 

It is clear from Table [4] that the Frob and Median split selection strategies are 
much better than the others for computing Frobenius numbers, and that Frob is a 
bit better than Median. We also see that applying the bound technique to the best 
split selection strategy improves performance by a factor of between two and three. 
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*: This time has been included despite using more than 512 MB of memory. 
Table 3. Empirical comparison of programs for irreducible decomposition. 



strategy 
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Table 4. Empirical evaluation of the bound technique. 
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